Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



What We Can Learn from Observations of Cooling Neutron 

Stars 



m 
o 
o 



> 

(N 
O 

a^ 
o 

o 

Oh; 

6 ■ 



D. G. Yakovlev^ and P. HaenseP 

^ loffe Physical Technical Institute, Politekhnicheskaya 26, 194021 St .-Petersburg, Russia, yak@astro.ioffe.ru 
Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland, haensel@camk.edu.pl 



Abstract. A generic toy model of a cooling neutron star (NS) is used to analyze cooling of NSs with nucleon 
and exotic compositions of the cores. The model contains the parameters which specify the levels of slow and fast 
neutrino emission as well as the lower and upper densities of the layer where the slow emission transforms into 
the fast one. The prospects to constrain these parameters from the present and future observations of isolated 
middle-aged NSs are discussed. 
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1. Introduction 



Neutron stars (NSs) are compact stellar objects which 
contain the matter of supranuclear density in their cores. 
The equation of state (EOS) of this matter cannot be cal- 
culated unambiguously (e.g., Lattimer & Prakash '2OUI5. 
Instead, there are many theoretical models which give a 
wide scatter of EOSs, from soft to stiff ones, with stan- 
dard nucleon/hyperon or exotic compositions of matter. 
The nucleon matter consists mainly of neutrons (n) with 
an admixture of protons and electrons (p and e) and possi- 
bly muons. Some model EOSs predict also the appearance 
of hyperons. The exotic matter may contain pion conden- 
sates, kaon condensates, or quarks (or a mixture of these 
components). The indicated models are still almost not 
constrained by the observations of NSs. 

In this paper we discuss the ability to constrain the 
EOS in the NS cores by confronting the observations of 
thermal emission from isolated NSs with the theory of NS 
cooling. In the last three decades the theory has been com- 
pared with the observations by many authors (e.g., Page 
I1998a|l but the problem is complicated and the main fea- 
tures are still unknown (too many factors are involved, 
such as neutrino emission mechanisms, superfluidity of 
baryon components of matter; see, e.g., Yakovlev et al. 
I2UU2|I . Thus, at this stage it is sufficient to be general 
rather than accurate and use a simple toy model of cool- 
ing NSs. It will allow us to formulate the current status of 
the problem (Sects. 4 and 5) without complicated simula- 
tions. 
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2. A toy model of a cooling NS 

Let us formulate our toy model. 

A middle-aged NS [t ^ 10^ yr) cools mainly via 
neutrino emission from its core (the region of densities 
p ^ 1.5 X 10"'^'* g cm~^). We assume that the core can be 
subdivided into three zones: the outer zone, p < ps] the 
transition zone, Ps l£ p < Pi] and the inner zone, p > pf. 
If the central stellar density pc < Ps, two last zones are 
absent. 

In the outer zone we assume a .slow neutrino emission 
while in the inner zone we assume a fast emission with the 
neutrino emissivity Q^, (erg cm~'^): 

™(P < Ps) = QsT|, Q'r'ip > pf) = QiTl (1) 

Here, Tg is the internal stellar temperature T expressed in 
10^ K, while Qs and Qi are constants. For simplicity, we 
use the linear interpolation in p between Q^J""" and Q^J"^^ 
in the transition zone. This approximation, although over- 
simplified, seems to be sufficient for exploring the main 
features of the transition zone - its position and thickness 
(which are presently unknown). 

The proposed generic description of Q^, covers a num- 
ber of physical model EOSs of nucleon and exotic supranu- 
clear matter with different leading neutrino processes col- 
lected in Tables 1 and 2. In these tables, N is a nucleon 
(n or p); v and i? are neutrino and antineutrino; q is a 
quasinucleon (mixed n and p states); u and d are quarks. 

For instance, Qs can describe modified Urea (Murca) 
process in nonsuperfluid nucleon matter, or weaker NN- 
bremsstrahlung (e.g., nn-bremsstrahlung if Murca is sup- 
pressed by a strong proton superfluidity as considered by 
Kaminker et al. I2UU11 12002(1 . The factor Qf can describe 
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Table 1. Main processes of slow neutrino emission in nu- 
cleon matter; Murca and bremsstralilmig (brems) 



Process 






Qs, erg cm s ^ 


Murca 


nN pNci? pNe - 


-+ nNz^ 


10^" - 3 X 10^^ 


Brems. 


NN ^ NNuiy 







the processes of fast neutrino emission: a powerful direct 
Urea (Durca) process in nucleon matter (Lattimer et al. 
I1991|l or weaker (but nevertheless strong) Durca-like pro- 
cesses in exotic phases of matter (pion condensed, kaon 
condensed, or quark matter) as reviewed, e.g., by Pethick 
H1992|l . The neutrino emission from hyperon matter is 
qualitatively the same as from nucleon matter. The bot- 
tom line of Table 2 refers to nonsuperfluid quark matter 
in NS cores. 

The transition zone mimics a onset of the fast neu- 
trino emission with growing p. In nonsuperfluid matter, 
the lower density ps is a threshold density of the fast 
emission; the threshold is usually sharp, i.e., pf ~ ps- 
In realistic models of superfluid matter, the fast emission 
turns on gradually, and the transition zone may be broader 
(e.g., Yakovlev et al. I2()()2(l . The broadening is caused by 
superfluidity provided its strength is high at the formal 
(nonsuperfluid) fast emission threshold and decreases at 
higher p. Then the superfluidity strongly suppresses the 
fast emission at the formal threshold and starts to open it 
at some higher p — ps, where the superfluid suppression 
ceases to be very strong. It fully opens the fast emission 
at still higher p = pf, where the superfluid suppression is 
almost completely removed. 

To be speciflc, we adopt a density profile within the 
NS in the form: p{r) = Pc[l - (r/R)^], where R is the NS 
radius. According to Lattimer & Prakash H2()()l|l . this is 
a reasonable approximation for many realistic NS models. 
Then the NS (gravitational) mass is M = 87rpc^'^/15. 

To follow the NS cooling we solve the equation of ther- 
mal balance in the approximation of isothermal interior 
(e.g.. Glen & Sutherland . 1980)1 : 

C{T,)^^-L^{T,)-L^{T,), (2) 

where Tg is the effective surface temperature, Ti{t) = 
T{r, t) e* is the redshifted internal temperature which is 
constant throughout the isothermal interior (p ^ pb ~ 
IqW g with account for the effects of General 

Relativity; T(r, t) is the local internal temperature of mat- 
ter, and $(r) is the metric function (describing gravi- 
tational redshift). Furthermore, C is the total NS heat 
capacity, L"^ = AiraT^ R"^ (1 — r^/R) is the photon sur- 
face luminosity as detected by a distant observer (rg = 
2GM/c^ being the gravitational radius), and is the 
redshifted neutrino luminosity. We have 

C{T,) ^ [dV c{T), L^{T,) ^ [dV Q,{T) e^*, (3) 



where c{T) is the heat capacity per unit volume, dV^ is an 
element of proper volume; integration is done over the NS 
interior. In addition, we introduce Tg°° = Tg y^l — r^/ R, 
the effective surface temperature detected by a distant 
observer. 

Taking the simplicity of the toy model, we assume 
dV — Airr^ dr (flat space) and e* = a/1 — rg/i? (constant 
redshift) in the NS interior. With the above assumptions 
on Qi, [T] , the function (T) is calculated in an analytic 
form. 

We employ c(T) = cq (T) /c , where cq (T) is the heat 
capacity of degenerate baryon matter composed of one 
particle species of number density n = p/m^ (ttin be- 
ing the bare nucleon mass) and effective mass mjj = 
0.7mN; /c is the parameter introduced to absorb the 
drawbacks of our toy model and to account for the ef- 
fects of other particles and superfluidity. The total NS 
heat capacity is then also evaluated in an analytic form: 
C = 47r X 0.2293 i?^Cc(T), where Cc{T) is the specific heat 
at the stellar center. In the NS models with Qs > 10^*^ 
(appropriate for nonsuperfluid nucleons, see above) we set 
fc = 1.25 adding thus 25% contribution of the heat capac- 
ity of protons fPaefe ll99^ . In the models with Qs < 10^° 
(appropriate to strongly superfluid protons at pc < Ps) we 
set fc = 1.0 at Pc < Ps and fc = 1.25 at pc > Ps- In fact, 
the value of fc weakly affects the cooling curves at the 
neutrino cooling stage and such variations of fc might be 
neglected: they do not change our principal results. 

We adopt the formula of Potekhin et al. H1997|l to re- 
late the internal and the surface temperatures (for the 
standard NS heat blanketing envelopes made of iron) . The 
surface gravity in this expression has been calculated as 
g = GM e""^ / R^ , i.e, including the General Relativity ef- 
fects. 

Thus we obtain a toy model of cooling NSs in a closed 
form. The model contains five parameters: Qs, Qi, Ps, Pi, 
and fc- Similar models have been used, e.g., by Lattimer 
et al. ipMjl and Page H1998al I1998b|l who have studied 
the sensitivity of the NS cooling to the efficiency of fast 
neutrino emission (i.e., to variations of Q{). They imple- 
mented the generic model of Q[f*'*(T) into exact cooling 
codes. Our toy model is much simpler but it enables us to 
make a general view of the problem (by varying addition- 
ally Qs, Ps, and pt) and mimic thus a number of physical 
phenomena without complicated computation. 

3. Cooling models 

For an adopted EOS of dense matter (for a set of five 
parameters, in our case) we can construct a sequence of 
NS models with different pc (i.e., different M). We can 
mimic the stiffness of EOSs by choosing different mass- 
radius relations. For simplicity, we take i? = 12 km for 
all models (the approximation of constant R may hold in 
a wide range of M for a number of EOSs, see Lattimer 
& Prakash 1200 1(1 . The main results will be the same for 
more realistic mass-radius relations. We will vary pc from 
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Table 2. Leading processes of fast neutrino emission in 
nucleon matter and three models of exotic matter 



6.5 



Model 


Process 






Qf, erg 


cm ^ s ^ 


Nucleon matter 


n ^ 




pe 


nv 


1026 




Pion condensate 


q 


qev 


qe 


qiy 


1023 


- 102« 


Kaon condensate 




qev 


qe ~* 


qiy 


1023 


-10^" 


Quark matter 


d ^ 


nev 


ue — > 


Au 


1023 


-10"* 



7 X 10"^'* g cm to 1.4 X 10^^ g cm ^ varying thus M from 
1.02 Mq to 2.04 Mq. 

The calculations predict three types of coohng NSs: [1) 
low-mass NSs, M < Ms (Mg corresponds to the central 
density pc = Ps)] [2) medium-mass NSs, < M < Mi 
(Mf corresponds to pc ~ Pf); and (5) high-mass NSs, M > 
M{. The same conclusions have been made by Kaminker et 
al. ( 2G02| with regard to the models of NSs with nucleon 
cores. 

Low-mass NSs are slowly cooling objects since their 
neutrino emission is slow: — Qf""". Their cooling 
curves are almost insensitive to M due to the obvious 
reason: both, L,y and C, increase with increasing M but 
this increase is compensated in the ratio L^/C which de- 
termines the cooling rate at the neutrino-cooling stage 
(L^ > L~f). 

High-mass NSs cool mainly via fast neutrino emission, 
Qy — Q^f*'*! from the inner zones. Their cooling curves are 
also not too sensitive to M . These rapidly cooling middle- 
aged stars are noticeably colder than the slowly cooling 
ones. 

Medium-mass NSs show the cooling intermediate be- 
tween the slow and fast ones. While increasing M from 
Ms to Mf and higher, we get a sequence of cooling curves 
which realize the transition from the slow to the fast cool- 
ing regimes. 

The thick and thin solid lines in Fig. ^ display typical 
cooling curves of low-mass NSs (pc = 8 x 10^"' g cm~^, 
M — 1.16 Mq) with two levels of the slow neutrino emis- 
sion, Qs- The long-dashed lines show cooling curves for 
high-mass NSs (pc = 1-4 x g cm'^, M = 2.04 Mq, 
Ps = 8 X 10"^^ g cm^'^, pf — 10^^ g cm^'^) with three levels 
of the fast neutrino emission, Qf = 10^^, 10^^ and 10^'', 
which roughly correspond to NSs with kaon condensates, 
pion condensates, and nucleon matter with Durca process 
in NS cores. NSs with hyperon cores are expected to cool 
at about the same rate as NSs with nucleon cores. One 
can see great difference of cooling scenarios at various Qs 
and Qf. 

The thick solid curve (Qs = 10^^) is very close to those 
obtained (e.g., Page ll998al I1998bl Kaminker et al. 2002) 
for low-mass non-superfluid NSs with the nucleon cores 
using exact cooling codes. It may be regarded as the basic 
slow- cooling curve. 

For comparison, we present three cooling curves calcu- 
lated with an exact cooling code for a 1.4 Mq NS (the nu- 
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6 7 



Fig. 1. Observational limits on surface temperatures of 
eight NSs versus toy-model cooling curves of low-mass 
NSs (solid lines) with two slow neutrino emission levels 
and high-mass NSs (long dashes) with three fast neu- 
trino emission levels. Dashed-and-dot line: slow cooling 
of a nonsuperfluid 1.4 Mq NS calculated with the exact 
cooling code; dotted line: the same but assuming triplet- 
state neutron pairing in the core, Takatsuka ( 1972) model; 
short-dashes: the same but with strong proton pairing and 
nonsuperfluid neutrons. 



cleon core with forbidden Durca process, pc = 1.22 x 10^^ 
g cm~'^, i?=11.65 km, EOS B in notations of Kaminker et 
al. l2()()2|l . The dot-and-dashed curve refers to a nonsuper- 
fluid star and agrees with the basic slow-cooling curve. 
The short-dashed curve is for a NS with strong proton 
superfluidity (which switches off the proton heat capacity 
and the neutrino reactions involving protons); it agrees 
with the upper solid curve. By varying the heat capacity 
parameter /c of our toy model, we could get even bet- 
ter agreement with the two exact cooling curves indicated 
above, ki t ^ 100 yr the exact code accurately describes 
the NS thermal relaxation which is not reproduced by the 
toy model. At these t, the exact cooling curves noticeably 
deviate from the toy- model ones. The third exact curve is 
explained later. 

Figure 13 shows the transition from slow to fast cooling 
with increasing M for the same Qs and Qf , ps = 8 x 10^^ 
g cm~^ (Ms = 1.16 Mq), and two relative density widths 
of the transition zone Ap/ps = (pf — Ps)/Ps = 0.01 or 0.1. 
The first value of Ap/ps corresponds to a sharp threshold 
of the fast emission, while the second value is appropriate 
to a smooth threshold. The NS age is fixed, t = 25000 
yrs (the age of the Vela pulsar, Lyne et al. I1996|l . The 
main features of Fig. |21 would not change if we varied ps 
from 8 X lO" to 1.2 x lO^^ g cm-^ (Ms from 1.16 Mq to 
1.75 Mq). If the contrast Qi/Qs between the fast and slow 
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0.8 1 1.2 1.4 
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Fig. 2. Transition from slow to fast cooling with increas- 
ing NS mass over the threshold value A/g for the same val- 
ues of Qs and Q{ as in Fig. ^ and two transition density 
intervals, Ap/ps = 0.01 or 0.1, at a fixed NS age t = 25000 
yrs. Shaded strip: observational limits for the Vela pulsar. 

emissivities is not too high {Qf/Qs ^ 10'^), the transition 
is rather smooth even for a very narrow density width Ap. 
For higher contrasts, the transition is smooth only if the 
density width is not too small, Ap/ ps <; 0.1. 

4. Confronting theory with observations 

Although the toy model is oversimplified, it reproduces the 
main features of exact calculations and can be confronted 
with the observations of thermal emission from cooling 
NSs. 

The observational basis is shown in Fig. ^ It displays 
the observational values of for six middle-aged iso- 
lated NSs, the same as in Yakovlev et al. ij^UU^ exclud- 
ing RX J1856-3754 and RX J0002-I-62. The latter two 
sources are most interesting but the interpretation of the 
observed spectra and the extraction of the surface ther- 
mal radiation is complicated (as discussed, e.g., in Pavlov 
et al. 120021 Walter & Lattimer |2)02 , and in references 
therein with regard to RX J1856-3754; as commented by 
Pavlov 2002 with regard to RX J0002-f62). The objects 
presented in Fig. dare: RX J0822-43, IE 1207-52, Vela, 
PSR 0656-hl4, Geminga, and PSR 1055-52. The data on 
Tg°° and t for these NSs are taken from the sources cited 
in Kaminker et al. 1)2002(1 . We display also the upper limit 
of for the Crab pulsar (Tennant et al. IM?T)l and for 
PSR J0205-H6449 in the supernova remnant 3C 58 (Slane 
et al. 120021. 

The comparison of the data with the cooling curves 
leads to the following generic conclusions. 

(1) The hottest observed NS, RX J0822-43, lies well 
above the basic slow- cooling curve but is compatible with 



the models of low-mass NSs with Qg 3 x 10^^, i.e., with 
the reduced slow neutrino emission. These conclusions 
have been made, e.g., by Kaminker et al. 1)20011 12002)1 . 
Alternatively, one can assume either additional reheat- 
ing mechanisms in the stellar interiors (which will com- 
plicate the theory) or the presence of accreted NS en- 
velopes (which increase the electron thermal conductivity 
and make the surface layers hotter, at the neutrino cooling 
stage, see, e.g., Potekhin et al. llH97l and Page ll99T^ . 

(2) The coldest observed NSs, particularly, the Vela 
and Geminga pulsars, lie well below the basic slow- cooling 
curve. They require enhanced neutrino emission in their 
cores, with Qf ^ 10^^, but we cannot pinpoint the na- 
ture of this enhancement. This conclusion has been made 
by several authors (e.g., Page ll99Halll998b| . Even a weak 
enhancement produced by kaon condensates would be con- 
sistent with the data, but any stronger enhancement pro- 
duced by pion condensates or nucleon Durca is also al- 
lowed. Thus, future search for colder NSs would be cru- 
cial but it may take time because of difficulties in detecting 
faint sources. 

(3) A rather uniform scatter of observational points be- 
tween the hottest and coldest NSs indicates a sufficiently 
smooth transition from slow to fast cooling, i.e., the ex- 
istence of representative class of medium-mass NSs. This 
circumstance has been mentioned, e.g., by Kaminker et 
al. pMl nna^ with regard to the coohng of NSs with 
nucleon cores. 

(4) The threshold density of fast neutrino emission, ps , 
can be placed anywhere in the interval from ^ 8 x lO^"' g 
cm"*^ to ^ 1.2 X 10"'^^ g cm~^. For any ps from this interval 
we can to build a sequence of models of medium-mass NSs. 
Tuning their masses, we could explain the data. Adopting 
various ps and pf we will attribute different masses to the 
same sources, as seen from Fig. 2 with the Vela pulsar as 
an example. Similar conclusions were made by Kaminker 
et al. 1)20011 |2D02) for coohng NSs with nucleon cores. 
Their transition layer was produced by the weakening of 
the proton pairing with increasing p and the associated 
broadening of the Durca threshold. 

(5) For a not too high contrast of slow and fast neu- 
trino emissivities (Qf/Qs ^ 10'^) we can build a repre- 
sentative class of medium-mass NSs even with a narrow 
density width of the transition zone, Ap/ ps <^ 1. For 
sharper contrasts, the density width has to be sufficiently 
wide, Ap/ Ps ^ 0.1, as follows from the cooling simula- 
tions of superfluid NSs with nucleon cores (Kaminker et 
al. 120011 12002)1 . One needs a wide transition zone to ob- 
tain a representative class of medium-mass NSs with su- 
perfluid nucleon or pion-condensed cores. This is possible 
in the presence of a strong superfluidity in the vicinity of 
the fast-emission threshold. 

5. Discussion and conc)usions 

Thus, our generic cooling analysis (5-parameter physics 
input leading to the families of cooling curves of NSs with 
different M) is too flexible to fix the nature of the fast 
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neutrino emission and its density threshold. Consequently, 
the data can be explained by a number of physical EOSs 
of dense matter. Notice that the cooling curves are more 
sensitive to the composition of matter than to the stiffness 
of the EOS (although the composition and stiffness are 
actually interrelated) . 

For instance, Kaminker et al. H2001[l2nn^ exploit the 
idea of nucleon matter in the NS core with the onset of the 
Durca process at high densities. They assume the presence 
of a strong proton superfluidity at not too high densities 
to suppress the Murca process and broaden the Durca 
threshold. The suppression of Murca allows them to re- 
duce the slow neutrino emission level Qs (from about 10^"'^ 
to about 3 X 10^^) and explain thus the hottest observed 
sources (Fig. The broadening of the Durca threshold 
ensures a representative class of medium-mass NSs to in- 
terpret cooler objects. 

In the phenomenological approach of Kaminker et al. 
l|^01 2002J the crucial ^Sq proton superfluidity is mod- 
eled by a specific density dependence of the proton critical 
temperature Tcp. Consider, for instance, a particular re- 
alization of this model applied by Yakovlev et al. ( 2002 ) 
to analyze the thermal state of PSR J0205+6449. For this 
specific model (model Ip with EOS A, in notations of the 
authors), the maximum of Tcp 7 x 10^ K) is reached 
at 5.75 X lO^"* g cm^'^. The critical temperature Tcp de- 
creases at higher densities, but it is still « 4.1 x 10^ K 
at the Durca threshold, pD = 7.85 x 10^"* g cm~'^, and 
then drops to zero at about 1.2 pD- Can Tcp be as high 
as ~ 4 X 10^ K at p = pd where the proton fraction is 
about 11%? In a recent paper Tsuruta et al. (.2002) ar- 
gue that this would contradict to the existing microscopic 
models of nucleon superfluidity, because the proton num- 
ber density Up is too large at p = po to allow for the 
""^Sq proton pairing. Taking the model of Yakovlev et al. 
H2002|l mentioned above we get Tip = 0.047 fm~'^ w 0.3 
at p = pd, i.e., noticeably smaller than uq = 0.16 fm~^, 
the nucleon number density in saturated nuclear matter. 
For the same number density of neutrons in neutron mat- 
ter, the critical temperature of the ^Sq neutron superflu- 
idity is ^ 1.3 X 10^° K if a realistic bare NN interaction 
is used (Lombardo & Schulze r2001|l . The medium effects 
(polarization effects, self-energy corrections) can decrease 
it to 7 X 10^ K or lower, depending on a particular model. 
In the case of the proton component of the neutron-star 
matter the np interaction can decrease the proton effec- 
tive mass, lowering further Tcp. However, this does not 
mean that the value Tcp(pd) ^ 7 x 10^ K, followed by a 
rapid drop of Tcp with increasing p, is ruled out by con- 
temporary microscopic theories. In contrast to the case of 
neutron matter at rin < tt-Oi reliable calculation of Tcp(pd) 
in the neutron-star matter, starting from a realistic NN in- 
teraction and including the medium effects, remains still a 
challenge for the many-body theory. Therefore, we think 
that the use of a "minimal model" (npe matter with Durca 
acting at p > Pd, combined with the appropriate nucleon 
superfluidity) remains a valid first-step approach. 



Notice that our general assumption on the neutrino 
emissivity (T, p) would be violated in the presence of 
('^P2) neutron superfluidity in the NS core with a den- 
sity dependent critical temperature Tcn(p) which has the 
maximum in the range from ^ 10^ to 2 x 10^ K (see, 
e.g., Kaminker et al. l2002|l . The neutrino emission due to 
Cooper pairing of neutrons will then be so strong that it 
will initiate a really fast cooling even at M < Ms, vio- 
lating the interpretation of relatively hot and old sources, 
first of all PSR 1055-52. As an example, by the dotted 
line in Fig. ^ we show the cooling of 1.4 Mq NS with the 
nucleon core and forbidden Durca process in the presence 
of a neutron superfluidity (model of Takatsuka ll972[ with 
maximum Tc of about 9 x 10* K at p « 5 x 10"^^ g cm~^). 
One can see that NSs with this superfluidity would be too 
cold to explain the data. 

Another physical model of cooling NSs was presented 
by Tsuruta et al. H2002|l . It is based on pion condensa- 
tion at supranuclear densities exploiting similar idea: slow 
cooling of low-mass NSs and faster cooling of massive 
NSs. Their cooling curve of low-mass NS (1.2 Mq, nu- 
cleon core, forbidden Durca process) agrees with our basic 
slow-cooling curve (after translating our curve, T^{t), to 
their format, L^{t)). Note that Tsuruta et al. I|2002|l mis- 
placed the positions of some observational data in their 
figure. Most important is RX J0822-43. They (as weU 
as we) take the data from Zavlin et al. H1999|l who give 
= (7.1 - 10.1) x 10^3 erg s'^ (IgL^ = 33.85 - 34.00) 
while Tsuruta et al. present IgT^ w 33.57 - 33.97. Thus, 
RX J0822-43 is sufficiently warm and cannot be explained 
with the basic slow-cooling model (Fig.^). Moreover, ac- 
cording to Tsuruta et al., they employ the model neu- 
tron superfiuidity of Takatsuka H1972|l . Therefore, their 
curve (if properly calculated) should resemble our dotted 
curve, and their model would then disagree with the num- 
ber of observational limits. To avoid this disagreement one 
can change the model of nucleon superfluidity. A natural 
model of rather strong proton superfluidity and weak neu- 
tron superfluidity at p ~ (3 — 8) x lO^'' g cm"'^ considered 
by Kaminker et al. H2001lOD02|l seems to be most suitable. 
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